Berry phase of magnons in textured ferromagnets 
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We study the energy spectrum of magnons in a ferromagnet with topologically nontrivial magne- 
tization profile. In the case of inhomogeneous magnetization corresponding to a metastable state of 
ferromagnet, the spin- wave equation of motion acquires a gauge potential leading to a Berry phase 
for the magnons propagating along a closed contour. The effect of magnetic anisotropy is crucial 
for the Berry phase: we show that the anisotropy suppresses its magnitude, which makes the Berry 
phase observable in some cases, similar to the Aharonov-Bohm effect for electrons. For example, 
it can be observed in the interference of spin waves propagating in mesoscopic rings. We discuss 
the effect of domain walls on the interference in ferromagnetic rings, and propose some experiments 
with a certain geometry of magnetization. We also show that the nonvanishing average topological 
field acts on the magnons like a uniform magnetic field on electrons. It leads to the quantization of 
the magnon spectrum in the topological field. 

PACS numbers: 75.45.+j; 75.30.Ds; 75.75.+a 



C3 



-a 
c 

o 
o 



> 
m 

o 

CO 

O 
in 
o 

-i— > 
a 



i 

C 

O 

o 



I. INTRODUCTION 

The Berry phase theorjiSi^ allowed to generalize the 
idea of Aharonov-Bohm effect^ on electrons in the elec- 
tromagnetic potential, to an analogous effect related to 
a gauge potential, which arises during the adiabatic mo- 
tion of a quantum system in a parametric space. Up 
to now a lot of efforts has been directed to understand 
better and to find an experimental confirmation for the 
Berry phase of electrons like, for example, in the case of 
electrons moving in a varying magnetization field of the 
inhomogeneous ferromagnet 

One of the most intriguing consequences of the Berry 
phase theory is a possibility of the Aharonov-Bohm-likc 
effect on electrically neutral particles or boson fields££ 
An example of the adiabatic phase for the polarized light 
has been investigated by Pancharatnam* and Berry^fi 
The other example is the Aharonov-Bohm effect for the 
excitonjii which is a bound state of an electron and a 
hole in semiconductors. 

Here we consider the effect of the gauge potential and 
Berry phase on the propagation of magnons in textured 
ferromagnets. Such quasiparticles are usually viewed as 
the elementary excitations of the ordered homogeneous 
state of a ferromagnet but they can be also used to clas- 
sify the excited states near a metastable inhomogeneous 
magnetic configuration. These magnons describe the dy- 
namics of weakly excited inhomogeneous ferromagnet. 

The dynamics of magnetization in nanomagnets is in 
the focus of recent activity^ because of the importance 
of this problem for magnetoelectronic applicationsii 3 ^ 
It includes the switching of magnetization by electric 
current, spin pumping, magnetization reversal in micro- 
scopic spin valves, etc. Usually, the magnons play a neg- 
ative role in the magnetization dynamics limiting the fre- 
quency of magnetic reversal, and also leading to the en- 
ergy dissipation. However, they can be probably used in 



the spin transport phenomena like the spin currents of 
magnetically polarized electrons. 

Here we study the energy spectrum of spin waves in 
ferromagnets with a static inhomogeneous magnetization 
profile, and we demonstrate a possibility of observation of 
the Berry phase in the interference experiments on spin 
waves in magnetic nanostructures. Recent results of the 
micromagnetic computer simulationi£ii£*i£ of such sys- 
tems demonstrate that the interference of spin waves can 
be really observed in magnetic nanorings with domain 
walls. 

The equation for spin wave excitations in a general 
case of arbitrary local frame depending on both coordi- 
nate and time, has been found long ago by Korenman 
et aii^ in the context of local-band theory of itinerant 
magnetismjia Here we use an idea of this method to relate 
the adiabatic space transformation to the Berry phase 
and to find corresponding properties of the spin waves in 
a topologically nontrivial inhomogeneous magnetic pro- 
file, which is a metastable state of the ferromagnet. We 
show that the magnetic anisotropy is a crucial element 
determining the possibility of observation of the Berry 
phase in real experiments. 



II. MODEL AND SPIN WAVE EQUATIONS 

We consider the model of a ferromagnet described by 
the Hamiltonian, which includes the exchange interac- 
tion, anisotropy, and the interaction with external mag- 
netic field . It has the following general form 



H 



dn 



^{n(r)} 



(1) 



where n(r) is the unit vector oriented along the magneti- 
zation M(r) at the point r, a is the constant of exchange 
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interaction, JF{n(r)} is a function determining the mag- 
netic anisotropy [correspondingly, it includes a certain 
number of tensors relating the components of vector n(r)] 
and the dependence on external field, and Mq is the mag- 
nitude of magnetization. 

Due to the condition n 2 (r) = 1, the model is con- 
strained and belongs to the class of nonlinear a models^ 
The stationary (saddle point) solutions for the magneti- 
zation vector n (r) describing metastable states of the 
ferromagnet, can be found by minimizing Hamiltonian 
(1) wi th the constraint n 2 (r) = 1. It was shown (see, e.g., 
Refs. [2l|22|) that such metastable states with inhomoge- 
neous magnetization profile are related to the topology of 
ferromagnetic ordering, and they can include skyrmions, 
magnetic vortices, and other topological objects. 

We will be interested in describing the dynamics of 
small deviations <5n(r) from a certain metastable profile 
no(r) with a nonuniform magnetization, n(r) = no(r) + 
<5n(r), |<5n(r) | <C 1. Correspondingly, we assume that the 
solution of a saddle-point equation describing the state 
no(r) is already known. 

We perform a local transformation 



n(r)=R(r)n(r), 



(2) 



using the orthogonal transformation matrix R(r). By 
definition, it determines the rotation of local frame in 
each point of the space, so that the magnetization in the 
local frame is oriented along the z axis, no = (0, 0, 1). 
Then we consider small deviations of magnetization s(r) 
from no. Since s(r) is small and vectors no are oriented 
along z, the vectors s(r) lie in the x-y plane. 

The transformation matrix in Eq. (2) is taken in a 
general form of orthogonal transformation 



R(r) 



e iip{r) J z e i9 (r) J y & i 



(3) 



where ip, 9, </> are the Eulcr angles determining an arbi- 
trary rotation of the coordinate frame, and J Xl J y , and 
J z are the generators of 3D rotations around x, y and z 
axes. 

Two rotation parameters (for definitcness, the angles 
9 and (f>) can be used to define the frame with the z 
axis along the vector no (r) . In the absence of anisotropy, 
the additional rotation to the angle ip is purely gauge 
transformation. However, in a general case of anisotropic 
system, this rotation allows to choose the local frame in 
correspondence with the orientation of anisotropy axes. 

The Hamiltonian of exchange interaction (the first 
term in Eq. (1)) in the rotated frame has the following 
form 



U - a 

ex ~ 2 



cPr 



dr.. 



Art 



(4) 



Transformation (3) and gauge potential (5) are 3x3 ma- 
trices acting on the magnetization vectors. The matrix 
Ai(r) can be also presented as 



A i (r) = iA?(r)J fl 



(6) 



where A^ (r) belongs to the adjoint representation of the 
rotation group. 

Using (3) and (6), we find the explicit dependence of 
the gauge potential on the Eulcr angles 

ax, \ . , 99 dcj> 

A,; (r) = sinw — sm9 cosw — — , 

on On 

d9 0(f) 
A\ (r) = cos ip 1- sin 9 sin ip — , (7) 

OTi OTi 

^(r) = - + costf-. 

The magnetic anisotropy described by the second term 
in the right hand part of (1) gives after transformation to 
the local frame a function ^{^(r)} with correspondingly 
transformed tensor fields. Here we do not restrict the 
general consideration of the problem by any specific form 
of the anisotropy but in the following we consider the 
most important examples of easy plane and easy axis 
anisotropy. 

The Landau-Lifshitz equations for the magnetization 
in the locally transformed frame are 



On 



SH 



SH 



dt - M ^ vxnu \5n x Vi S(Vfn T ) 



where e^x is the unit antisymmetric tensor, and 



d_ 

dr.; 



(8) 



(9) 



is the covariant derivative. The right-hand part of Eq. (8) 
vanishes for the magnetization profile corresponding to a 
metastable state. This is seen from the Landau-Lifshitz 
equation in the unrotated original frame. In the fol- 
lowing, we will use Eq. (8) for the small deviations of 
magnetization from the metastable state. Hence, we will 
consider in the right part of (8) only the terms linear in 
deviations. 

Using (1), (4) and (8) we find the equations for weak 
magnetic excitations near the metastable state (spin 
waves) 



ds x 
dt 



d 2 ' s v _ o A z — - (A z ) 2 i 



&4f 
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where the gauge field Ai(r) is defined by 



d 2 T 



Mq dn x dfiy 



{A\ys v + A*A\ S . 

j_ &jF_ 

+ M dfil & 



(10) 



( !; k R 1 



(5) 



dt 



drf 



2A*^L 
on 



(A 



z\2 r 
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dA z 2 
+ q t s v + (A\ ) Sx + A\ A i s y 

j_ d 2 T j_ 



(11) 



where c s = ja/Mo is the stiffness. 

Using (10) and (11) we can also present the equations 
for circular components of the spin wave, s± — s x ± is y , 



± i 



dt 



d 



=F i4 



On 
7 <9 2 .F 



V(r) 



7 a 2 jp 



2M dn 2 
s± + [-iy(r) - ic a .4f.4V 
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2M dn 2 
d 2 T 7 7 <9 2 .F 



M dn^dny 2M <9n 2 2M dhl 



(12) 



where V(r) and w(r) are, respectively, the effective po- 
tential and a mixing field acting on the spin wave: 



V(r)='f (4fr + (4fr 



I \ 8 



(4f f - (A\f 



(13) 



(14) 



Equations (12) for s + (r,t) and s_(r, t) are complex con- 
jugate to each other since they both describe the same 
spin wave with real components s x (r, t) and s y (r, t). 

We can see that V"(r) is an effective potential profile for 
the propagation of spin wave. Due to the terms w(r) and 
ic s AfA y in (12), the equations for circular components 
s + and s_ are coupled even in the absence of anisotropy. 
All these terms are of the second order in derivative of the 
rotation angle, and they are small in the adiabatic limit 
corresponding to a smooth variation of the magnetization 
vector n(r). 



III. SEMICLASSICAL APPROXIMATION 

Equations (10) and (11) can be solved in the semiclas- 
sical approximation. The condition of its applicability is 
a smooth variation of gauge potential 4f (r) and fields re- 
lated to the anisotropy, as well as the external magnetic 
field, at the wavelength of the spin wave, kL ^> 1, where 
k is the wavevector of the spin wave and L is the char- 
acteristic length of the variation of 4f(r) and ^ ? {n(r)} 
(more exactly, the minimum of the corresponding charac- 
teristic lengths). Note that the condition of applicability 
of the semiclassical approximation to solve the spin- wave 
equations, docs not require any smallness of the gauge 
potential itself. 

Starting from Eqs. (10) and (11), we look for a general 
semiclassical solution in the form 

s x (r, t) = a cos [£(r) — cut] + b sin [£(r) — u)t] , (15) 



s y (r, t) = d sin [£(r) — cut] + f cos [£(r) — cut] , (16) 

with arbitrary coefficients a, b, d, /, and a smooth func- 
tion £(r), so that we can neglect the second derivative 
of £(r) over coordinate r. Substituting (15) and (16) in 
(10) and (11), we can find four equation for the a, b, d, 
f coefficients. 

The solution (15), (16) describes the elliptic spin wave 
with an arbitrary choice of the axes x and y, and, gener- 
ally, with a varying in space orientation of the principal 
axes of the ellipse. We can simplify our consideration by 
choosing the angle ip(r) at each point of the space in ac- 
cordance with the orientation of the principal axes. The 
corresponding equation for ?/>(r) can be found from the 
condition of b = f = in Eqs. (15) and (16) 



Co A; Aa 



Mq dn x dhy 



0. 



(17) 



Using (17) and neglecting the terms with derivative of 
4f , which are small in the semiclassical approximation, 
we write the spin-wave equations (10) and (11) as 



ds x 
dt 



d ^-2Al d -^-{Alfs v+ {Ay? Sl 



+ M0 dh 2 y S 



V ' 



(18) 



dSy 

dt 



dr 2 



24^-(4?y.~ 



(At? 



_t_ a 2 ^ 

Mo ~dh 2 



(19) 



Note that by fixing the angle ip in Eq. (17), we are choos- 
ing the gauge, which defines completely the potential 4f . 
We do it in spirit of the usual fixing gauge in the WKB 
approximation. 

After substitution of (15) and (16) with b = f = into 
(18) and (19), we come to the following equation for the 
momentum fcj(r) = d£(r)/dri 

(u + 2c s At hf {c s [k 2 + (A!) 2 - (Af) 2 } 
+2 Px } {c s [k 2 + (At) 2 - (Af) 2 ] + 2 Py } - 0, (20) 



where 



7 d 2 T 



2M dfi 2 x y 



(21) 



are the anisotropy parameters. 

Equation (20) should be solved for fcj(r) func- 
tion of smooth inhomogeneous field 4f(r). This equa- 
tion does not constraint the orientation of k(r) but de- 
termines the magnitude of vector k(r) for each direction 
in the momentum space. Let us take vector k(r) along 
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an arbitrary direction, defined by a unity vector g. Then 
we can rewrite (20) as 

[w + 2c s kg^} 2 ~ {c s [kf + (At) 2 - (Af) 2 ] + 2p x } 

x{c s [k 2 + (Atf-(A!) 2 ]+2p y }=0, (22) 

and we come to the fourth-order algebraic equation for 
k(r). It can be solved numerically, and a resulting 
dependence of fcj(r) on the gauge field in the integral 
£(r) = J c fcj(r) dri leads to the Berry phase acquired by 
the spin wave propagating along the contour C. 

We can find the solution of Eq. (22) analytically in 
the limit of weak gauge potential \A\\ <C k, which cor- 
responds to the adiabatic variation of the magnetization 
direction n(r) and also the adiabatic rotation in space of 
the elliptic trajectory, \dip/dri\ <§; k. Then in the first 
order of Af we find 



>V2 



+ 9i9jA z j {r) (l+p 2 /u; 2 ) 



1/2 



-1/2 



(23) 



where p = \p x - p y \. 

Using Eq. (23) and taking the vector g along the tan- 
gent at each point of a closed contour C, we find the 
Berry phase 



1b(C) 



Afi^dn 



C (1+p 2 /lu 2 ) 



2^1/2 ■ 



(24) 



As follows from (24), the Berry phase 7b (C) in the 
anisotropic system acquires an additional factor k = 

1/2 

(l+p 2 /uj 2 ) depending on the magnetic anisotropy 
parameter p. 

The denominator in (24) has a simple geometrical in- 
terpretation. Indeed, the coefficients a and d in the semi- 
classical solution (15), (16) are the ellipse parameters, 
which are related to the anisotropy factor p 



= 1 



1/2 



P 

U> 



Correspondingly, we can relate the parameter n — 
(l + p 2 /u> 2 ) in Eq. (24) to the geometry parameters 
of the ellipse 



k = sin 29, 



(26) 



where = arctan(d/a). 

Using definition (7), the Berry phase can be finally 
presented as 



lB (C) = & [k - 1) 0) + « (cosfl - 1) £H dn 
Jc L o Ti ® r i 

(27) 

In this expression we extracted a term proportional to 
2nN , N e Z. This allows to avoid the multivaluedness 
of Berry phase in the absence of anisotropy when k = l. 23 



The first term in (27) is proportional to the total winding 
number of rotations associated with the angles tp and </>, 
whereas the second term is a spherical angle on 52 , which 
is the mapping space of the vector field n(r). The second 
term in (27) has a standard interpretation of the Berry 
phase as the magnetic flux penetrating the contour on 
1S2, when the field is created by monopole at the center 
of Berry sphere. Following this idea, one can interpret 
the first term in (27) as the flux created by the magnetic 
string along the z axis, penetrating through the mapping 
contour on the unit circlei^ 3 . In accordance with Eq. (27), 
this contribution to the Berry phase vanishes for isotropic 
magnetic systems, k = 1. The first term in (27) is the 
topological Berry phase (it depends only on the winding 
number) in contrast to the geometric Berry phase of the 
second term in (27) ~ 

As follows from (24), the effective gauge field for spin 
waves in the anisotropic system is A4 — kA%, and the 
corresponding topological field acting on the magnons 
can be calculated as the curvature of connection Ai (r) 



Bi = k e 



a4| 



-ijk 



Ok 

TT- A k ■ 



(28) 



Note that there is a contribution related to the variation 
in space of the anisotropy parameters [second term in 
Eq. (28)]. _ 

We consider now in more details the motion of ellip- 
tic spin wave in the adiabatic regime. The anisotropy 
suppresses one of the components s x or s y breaking the 
symmetry with respect to rotations around z axis. Cor- 
respondingly, there is no gauge invariance s + — > e ltp s+ 
and S- — > e^'^s- for the circular components, and the 
motion of magnetization in the spin wave is elliptical. In 
the adiabatic limit of \Af \ k, the solutions for s x and 
s y are given by Eqs. (15) and (16) with b = f = and 
the ratio of amplitudes (a/d). Thus, we could expect the 
local invariance to transformations preserving the value 
of s 2 + (d/a) 2 s 2 — const instead of simple rotations in 
x — y plane. 

Using the Fourier transformation of Eqs. (18) and (19) 

(25) 

for s x y we find the following equation for the elliptic 



components of spin wave, s± = s x ± i(d/a) s y 



2c s k>[A,j 



a 

2d 



= 0, (29) 



and the complex conjugate to (29), where c s = c s (1 + 
d 2 /a 2 )a/2d , and we determine the d/a from the condi- 
tion of vanishing of the second bracket in Eq. (29). This 
condition determines the ellipticity factor, and we find 
that it coincides with Eq. (25) in the limit of \Af\ <C k. 
Thus, we come to the following equation for the elliptic 
wave in the gauge field 

jw + 2c s hAl - 5 S [kf + (At) 2 ] - gj j 5+ = 0. (30) 
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This equation is not gauge invariant but in the adiabatic 
regime, neglecting the difference in small terms of the 
order of (Af ) 2 , we can present it as 



lo - c s (ki - nA*y 



pd 
2a 



S+ = 0. 



(31) 



Equation (31) contains a factor k < 1 before A\ and 
formally looks like the equation of motion of a particle 
moving in the reduced gauge field, which in turn leads to 
an effective suppression of the Berry phase. The calcu- 
lation of Berry phase using Eq. (31) with the gauge field 
suppressed by factor k leads us again to Eq. (24) . 

In the absence of anisotropy and in the adiabatic ap- 
proximation, the solution of spin wave equations has a 
simple form. The equations for circular components (12) 
are separated 



± 



d. 



..«») 



dt 



9 . 



,(o) 



and the corresponding solution is 



s+(r, t) ~ exp 



1 /2 

with ki(r) ~ giuj/cs +9igj Aj(r). The spin wave propa- 
gating along a closed contour C acquires the Berry phase 
of Eq. (24) with p = 0. 

Using Eqs. (7) we can present the topological field (28) 
in the absence of anisotropy as 



i / ki(r) dri — iojt 



(32) 



(33) 



Bi = — sin t 



eijk 



89 dcf> 
drj dr k 



(34) 



It does not depend on the angle ip, related to the choice 
of gauge like in the case of electromagnetism. 

By creating a certain metastable configuration of the 
magnetization no(r) in the ferromagnet, we simulate 
an effective gauge potential Ai(r), acting on the spin 
waves similar to the magnetic field in case of electrons. 
In particular, when the averaged in space topological 
field (28) is not zero, there arises the Landau quan- 
tization of the energy spectrum of magnons. In the 
absence of anisotropy, we find the quantized spectrum 
huj n = 2hc s |(B) I (n+1/2), where (...) means the average 
in space. 



IV. INTERFERENCE OF SPIN WAVES IN 
MESOSCOPIC RINGS 

Let us consider now the ring geometry of a ferromag- 
net with a topologically nontrivial metastable magneti- 
zation rio(r). It can be, for example, a magnetization 
vortex (Fig. la) or an even number of domain walls in 
one branch of the ring like presented in Fig. lb. Such 
a magnetization profile presents a metastable magnetic 
state. 




FIG. 1: Two rings with topologically nontrivial magnetization 
field. 



Let us consider first the case when there is no 
anisotropy. If kL <C 1 (adiabatic regime), the low- 
energy magnetic excitations of the metastable state are 
described by Eq. (32). Due to the presence of gauge po- 
tential A\ (r) , there is a phase shift of waves propagating 
from the point A, where the waves are excited, to the 
observation point B (see Fig. 1). The phase shift (Berry 
phase) equals to the integral § A\ (r) dri along the ring, 
and by using Stokes theorem can be calculated as the flux 
$ of topological field B defined in Eq. (34). It can be also 
presented as the spherical angle enclosed by the mapping 
of the ring to the circle at the unit sphere S% . This way 
we can find the phase shift of 2tt and irNd for Figs, la 
and lb, respectively, where an even Nd is the number of 
domain walls in the right arm of the ring. For example, 
in the case of Fig. lb with two domain walls in the right 
arm, there is no interference of spin waves excited in A 
and coming to the point B because the corresponding 
phase shift is 2ir. 

In the absence of anisotropy, the interference in the 
ring can be induced by rotating all magnetic moments 
from the plane to a certain angle a (the corresponding 
mapping is presented in Fig. 2b). The Berry phase associ- 
ated with the path along the ring will be smaller than 2ir. 
For a — 7r/6 the Berry phase turns out to be tt. It means 
that the experiment with interference of spin waves prop- 
agating from A to B through two different arms of the 
ring, would result in a complete suppression of the outgo- 
ing from B spin wave. Physically, it can be realized using 
the ring with very small easy-plane anisotropy, p/uj -C 1 
in a weak external magnetic field along z axis. 

A similar idea was recently proposed by Schiitz et aS^ 
for the radial orientation of magnetic moments under in- 
homogeneous magnetic field directed from some point at 
the axis of the ring. This magnetic field creates a "crown" 
of magnetic moments, and the corresponding mapping is 
similar to that presented in Fig. 2b. 

However, in the case of nonvanishing easy-plane 
anisotropy, there is no need to apply magnetic field 
to provide the interference of spin waves propagating 
in the geometry of Fig. lb. In this case the Berry 
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FIG. 2: Mapping of the ring to the n-space & (red contour) 
in the case of in-plane vortex magnetization shown in Fig. la 
(a) and for the same geometry with magnetization vector de- 
viating from the plane to the angle a (b). The Berry phase 
is 2-7T in case (a) and n for a — it/6 in case (b). 




FIG. 3: Spin waves on a wide ring (thin-wall cylinder made of 
easy-plane ferromagnetic ribbon) . The magnetization is along 
axis z. 



phase is given by Eq. (27) with p = ■yX/2M — const, 
and we obtain the difference in phases for two waves 
Ajb — itNd/(l + p 2 /uj 2 ) 1 / 2 . Thus, the interference of 
spin waves should be clearly seen for the two-arm ge- 
ometry with domain walls. The computer simulation 
experiments^eonfirm this expectation. 

Another possibility to observe the interference of spin 
waves can be presented by the geometry of a wide ring 
(thin-wall cylinder) like presented in Fig. 3. Assuming 
the easy-plane anisotropy of the ribbon, we obtain the 
ground state with a homogeneous magnetization along 
the axis of cylinder. The anisotropy axis is oriented radi- 
ally in each point of the cylinder, and the corresponding 



local frame is shown in the figure as x' y' z' . Due to the 
homogeneous magnetization, we get 9 = 0, and from (7) 
we obtain -4|(r) = d{ip + 0)/9rj. The components of 
anisotropy vector are u x = costp and u y = sin (p. Using 
Eq. (34) we find that the condition u y = reduces to 
ip + (j) = tp, and from (34) we obtain u x = 1. Correspond- 
ingly, the anisotropy parameter p = r y\/2Mo, and we 
obtain the Berry phase for the closed contour on the ring 
7b(C) = 2ttN c /(1 + p 2 /uj 2 ) 1/2 , were N c is the winding 
number of the contour C. 

The Berry phase of the spin wave propagating in mag- 
netic ring, plays the similar role as the phase of electron 
wavefunction in the Aharonov-Bohm effect with mag- 
netic flux penetrating through the ring. The role of mag- 
netic flux plays a string through the ring. 23 However, the 
flux created by the string does not depend on the size or 
shape of the magnetic ring. Correspondingly, the Berry 
phase associated with the string has the topological ori- 
gin, which makes it different from the Aharonov-Bohm 
effect induced by the magnetic-field flux thorough the 
conductive ring. 



SPECTRUM OF MAGNONS IN A RING 
WITH UNIAXIAL ANISOTROPY 



The other example is a ring with uniaxial anisotropy 
in a homogeneous magnetic field B along the axis z like 
presented in Fig. 4. Due to the anisotropy and exchange 
interaction, the magnetization along the ring is oriented 
like in Fig. 4, creating a certain angle 9 out of the ring 
plane. 

We take the anisotropy function T in the form 



(35) 



which corresponds to the uniaxial anisotropy along the 
ring. The local frame is chosen with the axis z along 
the magnetization at each point, and with the z-x plane 
tangential to the ring (parallel to the axis z). In this case 
we find T = — As 2 , cos 2 9/2. The angles determining the 
orientation of local frame are ip = and <fi = r/R. 

We can calculate the angle 9 using Eqs. (1) and (32) 
with vectors u and B along the axis z. Then, using the 
polar coordinates and the condition that n v and n z do 
not depend on the point along the ring, we find the energy 



/ an „ \n 
E = 2< i? U^f + ^ - M a Bn z 



(36) 



Substituting n v = sin (9, 



cos#, we calculate the 



angle 9 minimizing the energy (36) 

M B 



cos 9 = min < 1 , 



A - a/R 2 

for A > a/R 2 , and 9 = for A < a/R 2 . 



(37) 
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FIG. 4: Magnetization in a ring with uniaxial magnetic 
anisotropy in a homogeneous external magnetic field along 
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FIG. 5: Energy spectrum of magnons in the ring (Fig. 4) as 
a function of external magnetic field B, where K — A/2. 



The gauge potential in magnetic ring induces the en- 
ergy splitting of magnons propagating in the opposite 
directions^ Using Eqs. (7) we find A x = s'm9/R, A v = 
0, and A z — cos 9. It should be noted that the gauge 
potential is constant along the ring, so that there is no 
need to use the adiabatic approximation to determine the 
energy spectrum of magnons. 

Using Eqs. (10) and (11) after Fourier transformation 
over time t and coordinate r along the ring, we find 

UJS X = iC s [kl-Sy + 2ik n A Z S x + {A Z fSy + (A X ) 2 Sy] 

+ iA s y , 

LUSy = -ic s [k 2 s x - 2ik n A z s y + (A z ) 2 s x ] + ips x 

-iA s x , (38) 

where Ao = "fB. Here the momentum takes discrete val- 
ues k n = l-Knjh (n = 0, ±1, ...) to provide the periodicity 
of solution for the spin wave along the ring. We can use 
(38) to find the equation for elliptic components of the 
spin wave 



2 (lu + 2c s k n A z ) s + 

, c s (A*) 2 
H rip 



-c s {A z ) 2 

-c s {A z ) 2 

c s (A x ) 2 
V 



a ; 

Ao" 
VP 



1 



V 

7 2 -l 



(39) 



where rj = d/a is the ellipticity. The expression in second 
curved brackets vanishes for 



V 



c s k 2 n + c s (A z ) 2 + c s (A x ) 2 + &o 



"sk 2 



■ c s (A z ) 2 + A - p 



(40) 



Then, using (43) we find the energy spectrum of magnons 
in the ring 



27? 



A z 



VP 



4rj 



(A z 



^ s I Ax\2 _ 'If 

2 V [A ' 2 



2rj 



(41) 



The spectrum is shown in Fig. 5 for several first values of 
n as a function of the magnitude of field B. We take the 
parameters: A/2 — ttM 2 (it corresponds to the dipolar 
shape anisotropy of a magnetic cylinder), R — 75 nm, 
and XR 2 /a — 200. All the curves have a critical point B c 
corresponding to the magnetic field, for which the mag- 
netization M starts to deviate from the direction with 
9 = 0. In view of Eq. (37), B C M = A -a/R 2 . For n = 
the energy of magnons at this point is the soft mode with 
u> = 0. This mode corresponds to a uniform rotation of 
spins at each point of the ring toward a tangential to 
the ring direction. In the local frame it is the uniform 
rotation corresponding to the state with n = 0. 

For n and B — > we can find the spectrum in 
linear in B approximation 



w„(B) ~w„(0)+ 7 B 1 + 



1 



A\n\ (1 



,2)3/2 
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2c s nB 
\R 2 -a 



(42) 



where w n (0) = w_ n (0). It shows that the spectrum is 
degenerate at B = with respect to the sign of n. The 
splitting for i? ^ in linear approximation gives the 
curve for positive n below the one for the same negative 
in accordance with Fig. 5. 

At large magnetic field B, the systematics of levels 
should be changed. Namely, the lowest energy mode cor- 
responds to the uniform global deviation of orientation of 
all spins from the direction along the z axis. In the local 
frame, it corresponds to the mode with n = 1. Thus, it 
would be natural to label the modes with index m = n—1, 
so that the lowest in energy is the spin wave with m = 0. 
In the limit of B — > oo, each pair of modes ±m is degen- 
erate. It is clearly seen from Eq. (41) with r\ = 1. The 
splitting of these modes for B > B c demonstrates the ex- 
istence of topological Berry phase for the magnons on the 
ringS^ because the equilibrium state is the homogeneous 
magnetization, which leads to the vanishing of geometric 
Berry phase. 



VI. CONCLUSIONS 

We calculated the Berry phase associated with the 
propagation of magnons in inhomogeneous ferromagnets 



and mesoscopic structures with topologically nontrivial 
magnetization profile. We found that the most impor- 
tant effect is related to the magnetic anisotropy. Due 
to the anisotropy, the Berry phase for magnons is lower 
than a standard value of the spherical angle on the Berry 
sphere with a monopolc. Besides, an additional contri- 
bution to the Berry phase arises in anisotropic systems, 
which can be viewed as an effect of the gauge string pen- 
etrating through the mapping contour on the unit circle. 

Using these results, we demonstrated that the Berry 
phase can bc observable in interference experiments on 
magnetic rings, and controllable using the additional ho- 
mogeneous magnetic field. 
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